Direct current generation due to wave mixing in semiconductors 
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We describe a novel effect of the generation of direct current which may arise in semiconductors or 
' semiconductor microstructures due to a mixing of coherent electromagnetic radiations of commen- 

surate frequencies. The effect is, in essence, due to a nonparabolicity of the electron energy bands 
and is stronger in systems where this nonparabolicity is greater. We have made exact calculations in 
the framework of the Kane model, applicable to narrow gap semiconductors and the tight-binding 
model which we employ for a description of a semiconductor superlattice. 
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. The problem of emission and reception of electromagnetic radiation has attracted the attention of the scientific 
community for a long time. Good sources of coherent electromagnetic radiation, its receivers and detectors exist for 
the radio-frequency, microwave and optical ranges of the spectrum. Many of these devices are based on semiconductor 
technology. Nowadays, the terahertz range (0.1 - 1 THz, submillimetre wavelengths) is practically the last uncxploitcd 
. terra incognita. 

In 1970, Esaki and Tsu (!]] made a pioneering suggestion to use the semiconductor superlattice (SSL) for generation 
of Bloch oscillations of frequency ojb — eaEo/h, where Eg is a constant electric field applied along the axis of a SSL 
with a spatial period a. For typical SSLs and bias Eq oil - 10 kV/cm, the Bloch frequency belongs to the THz range 
PI . This can be used for the generation of THz radiation Q . The work of Esaki and Tsu fil stimulated enormous 
theoretical activity devoted to the interaction of a high-frequency electric field with a SSL Moreover, recent 

progress in the development of THz radiation sources and coupling techniques allow the systematic experimental 
' studies of many of the associated nonlinear effects ||. Amongst others, one of the most interesting suggestions was 
. that of Esaki and Tsu [[| and Romanov [Q to use the SSL as a new, artificial, nonlinear material for electromagnetic 
wave mixing and harmonic generation. The theory of wave mixing in SSLs, based on a solution of the Boltzmann 
equation, has been developed by Romanov and co-workers dfj]. 

Very recently, in cooperation with the Urbana group, we showed (|] that the effect inverse to the Bloch oscillations 
in SSL is possible, namely, an alternating field without constant bias can create Bloch oscillations in a single miniband 
^\ ' SSL. The appearance of induced Bloch oscillations means a spontaneous creation of constant voltage and corresponding 
0^ ■ direct current (DC) |l, that can be considered as a rectification of the THz field in SSLs. The rectification in || 
results from an interplay of dissipation (scattering of ballistic electrons with impurities and phonons) and generation 
of a self-consistent electric field along the SSL axis, which is a strongly nonlinear effect. Moreover, the generation of 
DC in § is a counterpart of chaos, which arises for the same ac field strength and frequency but for much weaker 
dissipation [ jTo| . 

Independently, Goychuk and Hanggi [Tl| ] suggested another scheme of quantum rectification using a wave mixing 
of an alternating electric field and its second harmonic in a single miniband SSL. The approach of pl[| is based on the 
theory of quantum ratchets (see |p^| and references therein) and therefore the necessary conditions for the appearance 
of DC include a dissipation (quantum noise) and an extended periodic system. 

In this work we find an even more general effect of DC generation, or, a rectification at wave mixing, which may 
arise in any semiconductor or semiconductor microstructure. The origin of this effect is related to neither dissipation 
nor quantum noise but rather to a nonparabolicity of the energy spectrum. This nonparabolicity is always present in 
any semiconductor or semiconductor structure Jl^Q ■ The values of the nonparabolical energy terms may take their 
largest values in semiconductors with wide bands or small effective masses. It seems that the size of the gaps is not 
important, however, semiconducors with wide gaps only have small concentrations of current carriers unless, of course, 
they are doped. Hence, the best candidates for the described effect are semiconductors having wide bands and narrow 
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gaps or doped semiconductors such as InSb and other analogous compounds. Indeed, wave-mixing has already been 
observed in such narrow-gap semiconductors p^-|l6|; its mechanism being mainly related to the nonparabolicity of 
the energy band @-|8|. 

We also consider another system, SSLs, which is now under intensive experimental investigation and the possibility 
of wave mixing there has also been suggested. The description of a SSL in the framework of the tight binding model 
allows exact solutions to be obtained. Moreover, since the nonparabolical energy term in a SSL takes smaller values 
than for normal semiconductors, the estimation of the DC generation effect for this system gives a lower bound to 
the value of the DC which may be observed in semiconductors. We show that the main features of the novel effect, 
the generation of DC due to the mixing of harmonics, may be well described and understood within the standard 
semiclassical Boltzmann equation approach . We demonstrate that the value of DC is strongly dependent on the 
product of a characteristic scattering time, r, and field frequency, f2. In the experimentally relevant case of weak field 
strength, the expression for DC can be represented in the form of the generalized Ohm's law with prefactor dependent 
on fir. The value of DC is maximal at fl ~ r _1 (namely, fir w 0.7) and it decreases quadratically for both small and 
large values of Sir. 

The novel effect we find is consistent with results obtained previously in a quantum approach |llf| for a special case of 
quantum noise. However, our predicted effect is much more general and has a much simpler physical picture associated 
with wave mixing in nonlinear media arising due to a nonparabolicity of the electron bands in semiconductors. In 
this new framework the results of [p|JTlJ|, which at first glance look very different, may be unified. The generality of 
the novel effect highlights new directions for various different applications and uses of semiconductor technology. 

Now, we demonstrate how the rectification effect arises due to a nonparabolicity of the spectrum of the semicon- 
ductors. Note that, although the dissipative property of a semiconductor or microstructure has an influence on the 
effect, it only plays a secondary role. To study the generation of a finite DC let us consider, for simplicity, a cubic 
semiconductor subjected to the electric mixing harmonic fields 



E(t) = Ex cos(f^) + E 2 cos(2fti 



(1) 



With the aid of an effective mass method jl3| the energy-momentum dispersion relation of a nondegenerate cubic 
semiconductor in the vicinity of the bottom of the conduction band may be represented as (see, for comparison, Ref. 

I) 
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pi, m is an effective mass at the bottom of the conduction band, r\ is a parameter of non- 



parobolicity. Following H, for simplicity, wc consider electron motion only along one direction, for example, along 
the x-direction and therefore, we may limit ourself to the one dimensional approximation. From now on the index 
x is suppressed. The electron's velocity is given by v = de/dp, the DC related to nonparabolicity is jd c — enrj(p 3 ), 
where n is the number of electrons per unit volume and the angled brackets (. . .) indicate the time averaging over a 
period of the electric field 2-7r/57. For pure ballistic electron motion without scattering by impurities or phonons, the 
electron dynamics is determined by the accelerating theorem p = eE(t). Combining this formula with Eqs. ([!]) and 
(0), we have 



j dc oc -enrj-—rElE 2 . 
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Next, we will consider two specific examples where the value of the DC can be expressed explicitly through the 
well-established parameters of the semiconductor energy bands. The energy-momentum dispersion relations, e(p), of 
a cubic semiconductor is usually obtained within the effective mass method (l3). For a SSL, the dependence e(p) 
was obtained within a tight- binding model H. In both cases, in the weak field limit when the excitation energy of 
electrons is small in comparison with the bandwidth, we may take into account only the first terms in the energy 
momentum relation (see, the dependence (||)). 

We start with the more transparent case of a SSL. Here the interplay of dissipation effects and nonparabolicity of 
the miniband may be described exactly. Consider the motion of an electron within a single miniband of the SSL with 
spatial period a and miniband width A. In the standard tight-binding approximation, we have 



A 



1 — cos 



(f) 



(4) 



2 



where p is the momentum of the electron along the SSL axis. For a weak electric field is the electrons oscillate in 
p— space near the center of the Brillouin zone \p\ <C irh/a or at the bottom of the miniband. Then, from Eq. (Q), we 
have (§) with r\ = -a 2 /6h 2 m and m = (2h 2 )/(Aa 2 ). The DC [Eq. (§)] takes the form j dc ~ f^-£i6, where £ ; = iff-, 
(i = 1,2). It is easy to see that the condition of weak nonparabolicity, \p\ <C nU/a, corresponds to <C 1. We turn 
to the detailed calculations of the DC taking into account collisions of electrons with impurities and phonons, i.e. 
dissipation effects. 

The electron transport properties in narrow miniband SSLs at temperatures above 40K are known pl[ to be well 
described by the scmiclassical Boltzmann equation with a constant relaxation time, r. Starting from a formal solution 
of the Boltzmann equation with constant relaxation time, Romanov and co-workers found the exact expression for 
a time-dependent current j(t) in a tight-binding lattice [Eq. (0)] subjected to an electric field with two frequencies 
U3\ and u)2 PJ22| (see also Appendix |A|). Taking oj\ = f2, 0J2 — 2Q and averaging j(t) over the period of the ac field, 
2tt/Q, we get for the rectified DC jdc = (j(t)) for t 3> t the following formula 

Ofi + 2/x 2 )a;cos(i/0) + sin(t/(/>) 

Jdc = JO 2^ 1 1 ~ 1 o..^2„.2 J M l(?l) J M2(?2)^l-2^(?l)^2+^(6), (5) 
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1 + (m + 2yL 2 ) 2 x 2 



where x — fir, jo — with a being a static, ohmic conductivity along the SSL axis, J M (£) is the Bessel function. 

The calculation of the DC with the use of Eq. (||) shows that the DC strongly depends on the product x = Sir. 
Fig. 1 illustrates the dependence of the DC on £1 and £2 for different values of x. For small x, the absolute value of 
the DC increases monotonically with an increase of both £1 and £2 (Fig. la). The increasing slope of the DC graph 
changes dramatically when the value of x increases. For instance, this slope can increase by almost five orders of 
magnitude with an increase in x of only one order (compare Figs, la and lb). However, for further change of a: from 
x = 0.2 (Fig. lb) to x — 1 (Fig. lc), the corresponding increase in DC slowsdown and becomes nonmonotonical. The 
DC reaches its maximal value for some optimal value of the relationship between £1 and £ 2 - So for the case x = 1, 
the DC is maximal when £1/^2 — 1.0 (see Fig. lc). 

When the electric field amplitudes are small, £1,2 <S 1, we can use the Bessel function approximation </„(£) ~ 
(£/2)™(l/n!) and obtain from (JsJ) the following analytic expression for the DC (see Appendix |b|) 

Jdc = - 1 - 5j0 Ax^+ X L 2 + l ^ 2COS<t> + O{e) - (6) 

Substituting usual Drude conductivity Odrudc = ne 2 T/m into jo, we get factor (enh) / '(am) , which is similar to 
the prefactor for jdc obtained within colisionless approximation. The equation (|^) agrees well with the asymptotic 
dependence of the DC on the field amplitudes, jd c oc £^£2 cos0 as obtained in reference |ll|] in the framework of a 
different approach. However, our result (^J) also gives the dependence of DC on the parameter x — Qt as well as 
indicate that the next significant contributions to the current appear only when the parameter, £, is of the order of 

r. 

Let us estimate the value of the predicted DC for the experimental conditions of THz field driven SSLs ||. For the 
miniband width A ~ 10 meV, superlattice period a ~ 10 nm, electron density n ~ 10 16 cm" 3 , sample area S ~ 10 
(/im) 2 , electric field of amplitude E ~ 1 kV/cm and frequency of several THz (£1 = £2 — 0.1) and the characteristic 
relaxation time of several picoseconds (x = Qt ~ 1), we get the DC Jdc = jdcS ~ 0.1 fiA. This value is in very 
good agreement with an estimate for the DC obtained by Goychuk and Hanggi fll| l from numerical calculations of 
the integral representations for DC within a different approach based on the ratchet theory. 

Note that DC density (|^) can be represented in the form of the generalized nonlinear Ohm's law incorporating only 
total electric field E(t) [Eq. ([!])]. Really, using the property (E 3 ) = 3/4(EfE 2 ) cos <f), we have from (||) the following 
expression for the DC 

*■—«*><*■*>. /M = fa<+ t + r (7 » 

The nontrivial prefactor f(x) increases as x 2 in the limit i< 1, reaches its maximal value f(x*) w 0.11 at x* w 0.71, 
and finally decreases as (4x 2 ) 1 for high frequencies (x = fir ^ 1). Note that in order to use the single miniband 
approximation and neglect the interminiband transitions, the electric field frequency, in units of the energy fiQ, should 
be less than the interminiband distance or, in other words, it should be of the order of or less than the miniband 
width HQ < A. The case HQ ~ A is the typical case in experiments || (see also our estimates above). Thus, the 
dimensionless electric field, E, (involved in the expression (^)) is of the order of eaE/A. To apply the Boltzmann 
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equation to the description of miniband transport in a quantum superlattice, the parameter eaE/A should always 
be small ||, i.e., in our notation, E <C 1. This remark shows that the DC generated in the SSL under the action of 
the high-frequency electric field E(t) is less than the corresponding current j — aE generated by a constant bias of 
the same strength E in the factor ~ 0.1E . The numerical value of this factor is only ~ 10~ 3 for typical experiments 
in SSLs, which implies that the semiclassical approach is valid and at the same time gives serious grounds for the 
observation of DC generation due to wave mixing. 

It is interesting that in another scheme of rectification as suggested in Ref. [Q only a single harmonic of a THz 
field is required. However, in this case there is actually an implicit wave mixing. Indeed, even in the simplest case 
of rectification, the self-consistent field generated in a SSL has a second harmonic of the fundamental frequency Q, 
so that the total electric field acting on the electrons has at least two components with frequencies f2 and 2f2 [p|,^0[. 
Consequently, the effect proposed in Ref. H may arise due to the wave mixing of the induced self-consistent field and 
the applied external THz field. This effect will be similar to the effect of DC generation arising at the wave mixing of 
two commensurate THz harmonics as described in Ref. PJll|. Thus, in spite of obvious differences, the approaches 
to rectification of THz radiation developed in Refs are not contradictory but rather complement each other. 

Now, let us consider the generation of DC at harmonic mixing in the III-V type semiconductors with narrow gaps, 
where the nonparobolicity is also strong. Mixing of mm-waves in n-doped InAs, InSb and GaAs semiconductors has 
been studied experimentally in ]l4|-|l6|] and discussed theoretically in |l7],|l8| a long time ago. However, we do not 
know of any reports devoted to the observation of DC generation at wave mixing in these semiconductors. This is 
probably because the experiments |l4|-[ll| were devoted only to the mixing of waves with similiar frequencies or third 
harmonic generation. 

Let us estimate the DC that could be generated in such semiconductors at harmonic mixing. Semiconductors such 
as InAs, InSb can be described @,H by the Kane four-band model @,||. Of course, for a proper description 
of the electron-hole dynamics we have to take into account all four bands, which is, in fact, a very tedious task. 
However, in a weak field assuming that, due to dominant donor doping, the electron concentration is larger than the 
hole concentration, we may take into account only one of the energy bands: the electron branch with the energy- 
momentum dispersion relation |19 |]l| in the form 



2p 2 n 1/2 



(8) 



where e g is the width of the gap. Note that GaAs cannot be described by the dependenc e (p|) , however in the weak 
field limit, the nonparabolicity, 77, is roughly twice as great as follows from the Kane model [|15( . Note also that, due to 
a diamond structure of the III-V type semiconductors, into a nonparabolicity coefficient the cubic invariants will also 
contribute, however, these contributions should not exceed those which follow from the Kane model. In the limiting 
case p/ yjme g -C 1, we have the dependence (||) with 77 = — l/(m 2 e g ) and the DC generated at wave mixing ([|) is of 
the order of 



The condition of weak field p/ '(mSg) 1 / 2 <C 1 takes the form Ei -C 1 (I = 1, 2). For the donor doped n-InSb with the 
gap width e g ss 0.2 eV, effective mass m sa 0.016m e , level of doping n ~ 10 16 cm -3 , sample area S ~ 1CP 4 cm 2 , field 
frequency ~ 10 11 sec -1 and field strength of ~ 100 V/cm, we get the DC idc = jdcS of the order of several mA. At 
the same time, the condition E ~ 0.1 <§; 1 guarantees the applicability of the weak field limit here. 

It is important to note that for gapless or narrow-gap semiconductors the current carriers may be excited by mi- 
crowave radiation. That is their concentration may be changed by the applied field. On one hand the nonparabolicity 
of the energy spectrum is larger in wide band semiconductors, while on the other hand the concentration of electrons 
which may be excited by microwave radiation may be larger in narrow gap or gapless semiconductors or semimetals. 
Therefore, the most optimal conditions for the described effect would be in semiconductors having narrow gaps and 
wide bands, like InSb and other analogous compounds. 

In this estimation of the DC in InSb we actually neglect a scattering of the electrons by impurities, defects and 
phonons; which is in itself an important process in semiconductors. In this respect, we should make two remarks. 
First of all, if the frequency f2 is of the same order as the inverse characteristic relaxation time r _1 , which is the case 

Otherwise, the dependence 



of experiments ]14|-|16|], the DC should be of the same order as estimated in Eq. (|9j) p 
of the DC will have a form similar to that obtained above for the case of a SSL, Eq. (J 



4 



Secondly, several other mechanisms of nonlinearity could also be responsible for wave mixing in semiconductors 
P^ , p7|jl8| . One of the most important of these is the so-called heating mechanism where the nonlinearity is related to 



the dependence of a relaxation constant on the field |16 18j]. This is mainly responsible for the mixing of mm- waves 
at low temperatures (~ 4K) |Q. On other hand, the experiments and theories [|l7 18 demonstrate that at 

temperatures ~ 80K the nonlinearity related to the nonparabolicity of the conduction band becomes more important. 
The estimates presented in |l8j show that the relative influence of these two different types of nonlinearity is dependent 
not only on the temperature but also on the type of wave mixing. For the mixing of waves with similiar frequencies 
(the case of experiments |l4j]), the nonparabolicity is the dominant nonlinearity for temperatures ~ 80K. However, 
for the mixing of waves with an integer ratio of frequencies, e.g. for a third harmonic generation, the influence of 
both types of nonlinearity is comparable |l6|,[l8| . In any case, the effect of DC generation arising due to the mixing of 
mm-waves should exist in narrow gap III-V and even some II- VI semiconductors. We hope that our results describing 
this new effect will attract the attention of experimentalists to this very intriguing issue. 

In summary, with the use of the semiclassical Boltzmann equation we have found the novel effect of DC gener- 
ation in semiconductors and semiconductor microstructures driven by two pure coherent electromagnetic waves of 
commensurate frequencies. The described effect originates from the nonlinearity associated with the nonparabolicity 
of the energy band and therefore, it is universal and may be observed in any semiconductors. Thus, our findings 
indicate that any semiconductor may be considered as a particular type of nonlinear medium serving as a generater 
of DC This effect is very different from the photogalvanic effect arising in semiconductors without a center of 

symmetry f24| . Finally, we propose that experiments be set up for the mixing of millimetre- and submillimetre- waves 
in semiconductors and semiconductor superlattices to observe the predicted generation of direct current. 

We thank David Campbell, Ethan Cannon, Laurence Eaves, Holger Grahn, Peter Kleinert, Eckehard Scholl, An- 
dreas Wacker for useful discussions on THz field driven semiconductor superlattices and Yuri Romanov for bringing 
our attention to refs. KNA thanks the Department of Physics, Loughborough University for hospitality and 

Krasnoyarsk Regional Science Foundation for partial support. 



APPENDIX A: 

In this Appendix we will derive Eq. (|^). Starting from the Boltzmann equation for a distribution function f{p,t) 

(Al) 



df | c£ df _ f-n 

dt dp t 



with an equilibrium distribution function f cq (p) and electric field E(t) — E\ cos(wii) -I- E 2 cos^i 
Romanov H have found the exact expression for time-dependent current 



j(t) =en v{p)f(p,t)dp 



in the form 



j(t)=j 



+00 



9 ^2 ^Mi (£0^2 (£2)^+1/1 (£1)^2+^2 (£2 

exp (-iviLu 2 t - w 2 {tjJ2t + 4>)) +c.c, 



Hl,p,2 = — OO v 1 ,l/ 2 = — OO 



1 + i([iiuji + [Iiu:i)t 



Orlov and 



(A2) 



(A3) 



where = ff^ (k = 1, 2) and jo — ^p. In this appendix, we first present simple derivation]] of Eq.(|A3|), and then 
will obtain the expression for rectified DC. 

Using periodicity of distribution function in p— space 



1 The simplest way is to solve the Boltzmann equation in the Fourier representation. For the case E(t) — Eocos(ut), this 
method is presented in clear form in: Bass F. G., Rubinshtein E. A., Fiz. Tverd. Tela, 19 (1977) 1379 [Sov. Phys. Solid State, 
19 (1977) 800]. 
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/(?,*)= E /«(<)exp(ix), 



the Boltzmann equation for ith harmonics could be represented as 



dfi iela dA ft- /, 

Jl = — 
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9i Tic 9i " r 
Assuming /(p, f = — oo) = f cq (p), the solution of (|Aq) is 
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(A4) 



(A5) 



(A6) 



This solution is valid for arbitrary A(t); in our case A(t) = —^jf sin(wii) — sin(w2i + 0). Using the Bessel functions 
expansion 



+00 



exp(±«£sin#) = E Jm(Q exp(±im#), 



(A7) 



m— — 00 



we can expand exp[zb 



ielaA(t) i 



in (AG) and get after integration: 



/i = /r q E E ^(6)^(6)^(6)^(6)- 



i(n 1 -k 1 )u 1 t i(n 2 -k 2 )uj 2 t 

e i(n 2 -k 2 



ki,k 2 n 1 ,n 2 



(A8) 



Now, introducing new summation indexes Vi = ki — ni (I = 1, 2) and substituting Eq. ( A8) into the expression for the 
current j = (/1 — we get the formula (A3). The Fourier transform of equilibrium distribution function 

determines the temperature dependence of conductivity. 

Finally, taking w\ = SI, W2 = 2f2 and averaging j(t) [Eq. (|A3|) j over period of ac field 2n/£l, we get for the rectified 
DC j dc = (j(t)) the formula (§). 



APPENDIX B: 

In this Appendix we derive Eq. @. We start with the expression for DC jdc, Eq. (JsJ) , which now is convenient to 
rewrite in the form 

+00 +00 

jdc = jo E E [^(/"i'M2,^) + 5(mi, H2,v)\ J m (6)^2(6)^1-2^(6)^2+^(6), (Bl) 

MI1M2 — — 00 1/ — — oc 

where 

(Mi +2/x 2 )xcos(t/(/)) sm(t^) 

^(Mi.M2,f) = -r— 7 — ^T' S(mi,M2,^J = r— 7 : , 2 2 (B2) 

1 + (Mi + 2 M2j 1 + (Mi + 2M2) £ 

Using the Bessel function property J n (£) ~ 6 f° r £ ^ 1; we wm consider only nonvanishing lowest order in 6,2 <C 1 



terms in (Bl). We start with terms of [mi = 0, M2 = 0, v = ±1]: 



[A(o,o, 1) + s(o,o, 1)] J (6)</o(6) J- 2 (6) + 
[A(o, 0,-1) + s(o, 0,-1)] j (6)Jo(6)^(6)^i(6) = 

2 sin(</>) J (6 ) Jo (6 ) J 2 (6 ) Ji (6 ) • (B3) 



Term with [0, —1, 1] plus term with [0, 1, —1]: 



[A(o, -1, 1) + s(o, -1, 1)] J (6) J_ 2 (£0 Jo(6) + 

[A(0, 1, -1) + 5(0, 1, -1)] Jo(6)^i (6)^2(6) Jo(6) = 
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2x cos d> sin < 



1 + 4x 2 1 + 4x 2 
Term with [1,0,1] plus term with [-1,0,-1]: 



[4(1,0, 1) + 5(1,0, l)]Ji(ti)M&)J-i(Si)Ji{t2) + 
[A(-l, 0,-1) + B(-l, 0, -1)] j-i(a) Jo(6) Ji(6) 



x cos ( 



1 + x 2 1 + x 2 
Term with [1,-1,1] plus term with [—1,1,-1]: 



■J? (6)^0(6) Ji (6). 



[4(1, -l, i) + -l, i)]Ji(Si)J-i{&)J-i(Si)MS2) + 

[a(-i, i, -l) + b(-i, i, -l)] Jb(6) : 



— x cos cj> 
1 + x 2 



Term with [2, 0, 1] plus term with [-2, 0,-1]: 



[4(2,0, 1) + 5(2,0, 1)] J 2 (£i) J (6) Jo(£i)Ji(&) + 
[A(-2,0,-l)+ J B(-2,0,-l)]J_ 2 (a)Jo(6)^o(eOJ-i(6) 



2x cos < 



1 + 4x 2 1 + 4x 2 
Term with [2, —1, 1] plus term with [—2, 1, —1]: 



J2^l)Jo(6)Jo(6)Jl(6)- 



[4(2,-1, 1) + B(2, -1, l)]M£i)J-ifa)Mti)Jifa) + 
[4(-2, 1, -1) + B(-2, 1, -1)] J-2(6)^i(6)Jb(6)Jb(6) 



(B4) 



(B5) 



(B6) 



(B7) 



-2am(4>)J2(Si)Ji(&)MZi)MS2 



(B8) 



All other terms in Eq. (Bl) cither compensate each other or are of order which is greater than £ 3 . Combining 
(p3|)-(p^), we get following formula for the DC 



1 + X 2 



1 



£i6cos( 



2 V 1 + 4x 2 1 + x 2 
where we have used the expansion 

J n (0 « K/2)"(l/n!), £<1. 



2 4x 4 + 5x 2 + 1 



cos </>, 



(B9) 



(B10) 



As follows from Eq. (B9), the absolute value of the DC as a function of x = VLt at fixed £i, £2 and <f> increases like 



x 3 for small x, then reaches maximum at x* » 1.3 {\jdc(% — x *)/jo\ ~ 0.16£ 2 £2 cos</>), and finally decreases as 3/(8x) 
for large x. 

To check these calculations and understand what is the order in £^2 for the next terms in expression for DC, we 
perform the symbolic computations for Eq. (Bl) with account of ( B10| ) using MAPLE. As example, we present here 
the form of DC obtained by the summation in (Bl) over intervals —2 < (/zi, /i 2 , v) < 2: 
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1 



sva((j)) 



589824 1 + 4ai 2 
1 sin(0) 



(1 cos(</>) x 
294912 l + 4x 2 

(1 cos(</>) a; 
4096 1 + 36 a; 2 

(1 cos(^) X 
768 1 + 16x 2 ~ 3072 1 + 16.x ; 

1 cos(</>) a; 1 sin(0) \ 



i 



24576 1 + 36 a; 2 
1 sin((/>) 



24576 

1 cos(0) X 



1 



sin(<^>) 

1536 1 + 4 a; 2 ^ 3072 l + 4a; 2 



6144 l + 9a: 2 18432 l + 9a; 2 



s5 



/ 1 cos(0) a; 1 sin((/>) 5 cos(</>) a; 1 sin(^>) 
\ 768 1 + x 2 + 768 1 + x 2 + 768 1 + 25 a; 2 ~ 768 1 + 25 a; 2 

/ 1 cos(0) a; 1 sin(^) 1 cos((f>) x 1 sin(</>) 



+ 



32 l + 9x 2 
1 sin(</>) 



96 1 + 9 a; 2 96 1 



5 cos(^) x 
1536 l + 36x 2 ~ 768 l + 25x 2 

/ 3 cos(</>) x 1 cos(</>) x 



: 2 96 1 
sin(</>) 



1 

768 1 + 25 a; 2 
1 cos(</>) a; \ 



16 1 + 9a; 2 8 1 + 16a; 2 16 1 + 4x 2 



1 cos(</>) a; 
384 1 + 16x 2 " 

/ 1 cos(</>) a: 




1 cos(^) a; 
256 l + 36x 2 

1 sin(^) \ 



1536 1 + 16 a; 2 
1 cos(^) X 



2 1 + 4x 2 2 1 



6 Kr 



As is evident from this result, next after (B9) contribution do DC is of order of £ 5 (i.e., of order of £{£2 with /c + Z = 5) 
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FIG. 1. The dependence of DC on amplitudes E\ and E2 [Eq. (5)]. The values E\ and E2 are expressed in units of hu>/(ea) 
and 2%w/(ea) 1 respectively, and DC is normalized to jo. The product Qt has the values: x = 0.01 (a), x = 0.2 (b), and x = 1.0 
(c). Phase = 0. 
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